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This  paper  develops  a self-consistent  theory  of  the  cyclotron 
maser  instability,  assuming  azlmuthally  symmetric  perturbations  about 
a slowly  rotating  hollow  electron  beam  propagating  parallel  to  a uniform 

A 

axial  magnetic  field  Bq£z>  The  stability  analysis  is  carried  out  within 

the  framework  of  the  linearized  Vlasov-Maxwell  equations.  It  is  assumed 

that  the  beam  is  thin  with  radial  thickness  (2a)  much  smaller  than  the 

2 2 

beam  radius  (Rn) , and  that  u>  /u>  <<1,  where  cu  and  u are  the  electron 
U pc  pc 

plasma  frequency  and  electron  cyclotron  frequency,  respectively,  in  a 
frame  of  reference  moving  with  the  beam  axial  velocity  cB^-  The  analysis 
is  carried  out  for  the  specific  choice  of  equilibrium  electron  distribu- 
tion function  in  which  all  electrons  have  the  same  value  of  canonical 
angular  momentum  and  the  same  value  of  energy  in  a frame  of  reference 
moving  with  axial  velocity  cB^.  Stability  properties  are  investigated 
including  the  important  Influence  of  finite  radial  geometry,  finite  beam 
temperature,  and  transverse  magnetic  perturbations  (6B^0).  It  is  shown 
that  instability  exists  for  a very  narrow  range  of  axial  wavenumbers 
satisfying  | k-B^u/c | <<1/Rq.  Detailed  stability  properties  are  calculated 
for  a variety  of  system  parameters. 
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I.  INTRODUCTION 

There  is  a growing  literature  on  the  equilibrium^  ^ and  stability^  ^ 

properties  of  intense  hollow  electron  beams,  including  the  recent  experimental 

19 

studies  of  basic  stability  properties,  plasma  confinement  schemes 

20  21-25 

such  as  Astron,  high-power  microwave  generation,  and  electron 

25 

ring  accelerators.  Three  of  the  basic  instabilities  that  occur 

5 6 19 

in  a hollow  electron  beam  are  the  diocotron  instability,  ’ * the 

8-11  21-24 

electron  cyclotron  maser  instability,  * and  the  negative-mass 

12—18  25 

instability.  * The  negative-mass  instability,  which  is  an  important 

instability  for  a fast  rotational  equilibrium,  is  associated  with 
azimuthal  charge  bunching.  On  the  other  hand,  the  electron  cyclotron 


maser 


8—11  21—24 

instability,  ’ which  is  a dominant  instability  for  a slow 


rotational  equilibrium,  is  present  even  for  azlmuthally  symmetric  pertur- 


bations with  9/30-0.  However,  the  physical  mechanism  for  the  cyclotron 


maser  instability  is  similar  to  that  for  the  negative-mass  instability. 


14-18 


with  radial  charge  clumping  dominating  for  £-0.  Here  l is  the  azimuthal 
harmonic  number.  The  free  energy  for  driving  the  instability  is 
associated  with  the  transverse  energy  spread  of  the  beam  electrons. 

This  paper  develops  a self-consistent  theory  of  the  cyclotron  maser 

instability  for  azlmuthally  symmetric  perturbations  about  an  intense 

hollow  beam  equilibrium.  The  analysis  is  carried  out  within  the  framework 

of  the  Vlasov-Maxwell  equations  for  an  infinitely  long  beam  propagating 

parallel  to  a uniform  magnetic  field  with  axial  velocity 

(Fig.  1).  The  positive  ions  form  an  imnobile  (m^-w*)  partially  neutralizing 

background.  In  addition,  we  assume  that  the  electron  beam  is  thin  [Eq.  (1)], 

i.e. , the  radial  thickness  (2a)  is  small  in  comparison  with  the  mean 

2 2 

radius  (Rn).  It  is  also  assumed  that  u <<<d  , where  u>  and  u are  the 
« p c p c 

electron  plasma  frequency  and  electron  cyclotron  frequency,  respectively. 


in  a frame  of  reference  moving  with  axial  velocity  V,^.  Equilibrium  and 
stability  properties  are  calculated  for  the  specific  choice  of  equilibrium 
electron  distribution  function  [Eq.  (8)], 

f0fu  _ . noRo 

fe(H’P0’Pz>  


2 SW&Gq-Vq)  . 


2lnnVo 

where  H is  the  energy,  P.  is  the  canonical  angular  momentum,  p is 

O Z 

2 

the  axial  momentum,  U=H-6^cpz~YQmc  /y^  is  an  effective  energy 
variable,  and  Uq , Rq,  Yg»  y^»  8^  and  Pq  are  constants.  Equilibrium 
properties  are  examined  in  Sec.  II.  One  of  the  important  features  of 
the  analysis  is  that  the  equilibrium  distribution  function  in  Eq.  (8) 
corresponds  to  a sharp-boundary  density  profile  [Eq.  (17)],  with  uniform 
axial  velocity  over  the  beam  cross  section  [Eq.  (18)],  and  nonzero  axial 
and  radial  beam  temperatures  [Eq.  (21)]. 

Stability  properties  are  investigated  in  Secs.  Ill  and  IV,  including 
the  Important  influence  of  (a)  finite  axial  and  radial  temperature, 
which  limits  the  unstable  range  of  k-values  to  a narrow  bandwidth 
(Appendix  C) , (b)  finite  radial  geometry  and  the  presence  of  a grounded 
conducting  wall  at  radius  r”Rc*  and  (c)  the  influence  of  transverse 
magnetic  perturbations  (6^0).  As  shown  in  Sec.  Ill,  the  most  unstable 
perturbation  is  the  TE  mode  with  axial  wavenumber  k-S^oi/c.  Here  u>  is 
the  eigenfrequency,  and  cg^EV^  is  the  axial  velocity  of  the  beam. 

The  stability  analysis  predicts  Instability  for  a very  narrow  range 
of  k values  corresponding  to  j kc— | <<c/Rq  (Appendix  C) . In  this  context, 
the  TE  mode  dispersion  relation  can  be  approximated  by  [Eq.  (71)] 


p (b_+b+)+(u>-u)0) 


(b-+V 


v(a)“koV 


“0’k0  VoV^Vb-VV 


1 - 


wy 


^““VWV  J 


for  kc*8v(o.  Here,  b *[b  ] , are  the  normalized  magnetic  wave 

D X ± 0)« , 

u u *2  2 2 2 
admittances  defined  in  Appendix  B,  v is  Budker's  parameter,  p "<*>q/c  -kg, 

2 1/2 

Bq-CYq-!)  /yq»  and  (aig,kg)  is  defined  in  Eq.  (66).  A detailed 


numerical  analysis  of  the  dispersion  relation  (71)  is  carried  out  in  Sec.  IV. 
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A A 

It  is  found  that  the  geometric  configuration  (b  +b+)  and  electron  density 
(v)  have  a strong  influence  of  stability  behavior.  For  low  beam  densities, 
instability  exists  only  for  (b_+b+)>0. 

21 

As  a comparison  with  previous  analyses,  it  is  instructive  to 
consider  the  limit  of  a very  tenuous  (v/yq-H)) , infinitesimally  thin 
(a/Rg-*0)  beam.  In  this  case,  Eq.  (71)  can  be  approximated  by  [Eq.  (69)] 


5 n2 

u>  .2  °0n  2v 

Jl(a0nW 

2 

'(u)-kVb) 

^(m-kv/ 

c 

RcJ2<“0„> 

(co-kVb-a)c/Yb) 

3(u)-kVb-a)c/Yb)2 

where  o.  is  the  nth  root  of  J, (a_  )-0  and  J„(x)  is  the  Bessel  function 
On  1 On  l 

of  the  first  kind  of  order  l.  The  general  features  of  Eq.  (69) 

are  qualitatively  similar  to  the  dispersion  relation  obtained  by 
21 

Friedman  et  al.  for  perturbations  about  a self-consistent  Vlasov 
0 2 

equilibrium,  fe*const.’<6(H1-YQmc  )6(pz-Y^m6^c)6(P0-PQ)  , where  all 

2 

electrons  have  the  same  perpendicular  energy  H^Yq111^  • However, 

Eq.  (69)  differs  in  detail  as  do  the  corresponding  stability 

results.  For  example,  the  final  term  on  the  right-hand  side  of 

Eq.  (69)  is  proportional  to  Bn(o>-kV,  ) 2/3,  instead  of  Bn(u^-k^c^) /2 

u b u 

21 

as  obtained  in  Ref.  21.  Moreover,  in  the  previous  study,  the 
system  can  be  unstable  for  a broad  range  of  axial  wavenumbers  whereas 
the  present  stability  analysis  predicts  a narrow  band  of  unstable  k- 
values . 
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II.  EQUILIBRIUM  THEORY 


A.  Equilibrium  Configuration  and  Basic  Assumptions 

The  equilibrium  configuration  is  illustrated  in  Fig.  1.  It 
consists  of  a nonneutral  hollow  electron  beam  that  is  infinite 
in  axial  extent  and  aligned  parallel  to  a uniform  applied  magnetic 

A 

field  The  mean  radius  of  the  electron  beam  is  denoted  by  R0 

and  a grounded  cylindrical  conducting  wall  is  located  at  radius  r-R  . 

c 

The  applied  magnetic  field  provides  radial  confinement  of  the  electrons, 
and  the  radial  thickness  of  the  electron  beam  is  denoted  by  2a. 

As  shown  in  Fig.  1,  we  Introduce  cylindrical  polar  coordinates  (r,0,z) 
with  z-axls  coinciding  with  the  axis  of  symmetry;  r is  the  radial 
distance  from  the  z— axis,  and  0 is  the  polar  angle  in  a plane  perpendicular 
to  the  z-axls.  The  electron  charge  is  partially  neutralized  by  a 
positive  ion  background. 

The  following  are  the  main  assumptions  pertaining  to  the  equilibrium 
configuration: 

(a)  Equilibrium  properties  are  independent  of  z (3/3z-0)  and 
azimuthally  symmetric  (3/30-0)  about  the  z-axis. 

(b)  The  mean  canonical  angular  momentum  of  the  electrons  is 
negative,  which  corresponds  to  a slow  rotational  equilibrium.4,6 
The  positive  ions  form  an  immobile  (m^**),  partially  neutralizing 
background . 

(c)  The  radial  thickness  of  the  electron  beam  is  much  smaller  than 
the  beam  radius,  l.e.. 


a«R 


0 * 


(1) 


(d)  It  is  further  assumed  that 


6 


v a , 2 ..  . 2 

^ V (y<T1)/yo  • 

2 2 

where  v»N  e /me  is  Budkerfs  parameter, 
e 


r 2 

N =2ir 

S V 


dr  r n (r) 
e 


is  the  number  of  electrons  per  unit  axial  length,  ng(r)  is  the  equilibrium 

electron  density,  c is  the  speed  of  light  in  vacuo,  -e  and  m are  the 

2 

charge  and  rest  mass,  respectively,  of  an  electron,  and  YQmc  is  the 

electron  energy  in  a frame  of  reference  moving  with  the  mean  axial 

velocity  V^z  of  the  electron  beam.  In  Eq.  (3),  and  denote  the 

inner  and  outer  radii,  respectively,  of  the  beam.  The  inequality 

in  Eq.  (2)  can  be  expressed  in  an  alternate  form.  For  a thin  beam, 

from  Eq.  (3),  the  average  electron  density  ng  is  approximately 

n=N  /(4irR_.a).  Combining  this  result  with  Eq.  (16),  Eq.  (2)  can  be 
U e u 

expressed  as 


(i>2/oj2«1  (4) 

P c 

2 1/2 

where  u “(4irarte  /Y-.m)  and  oj  meBn/ynmc  are  the  electron  plasma 
p u u c u u 

frequency  and  electron  cyclotron  frequency,  respectively,  in  a 
frame  of  reference  moving  with  axial  velocity  V^. 

(e)  Consistent  with  the  low-density  assumption  in  Eqs.  (2)  and 
(4),  we  neglect  the  Influence  of  the  small  equilibrium  self-electric 
field  E®(r)j£r  and  self-magnetic  field  B® (r)^g+Bz(r)^z  that  are  produced 
by  the  lack  of  equilibrium  charge  and  current  neutrality. 


B.  Hollow  Beam  Equilibrium  Properties 


Central  to  a description  of  steady-state  Vlasov  equilibria 
are  the  single-particle  constants  of  the  motion  in  the  equilibrium 


7 


field  configuration.  For  azimuthally  symmetric  equilibria  with 
3/30*O“3/3z,  there  are  three  single-particle  constants  of  the  motion. 
These  are  the  total  energy  H, 


„ 2 , 2 4 2 2.  1/2 

H-ymc  = (m  c +c  £ ) 


(5) 


the  canonical  angular  momentum  P , 

0 


p0=r [p0-(e/2c) rB0 ] » 


(6) 


and  the  axial  canonical  momentum  P , 

z 


Pz=P2  * 


(7) 


where  the  equilibrium-self-fields  have  been  neglected  in  comparison 
with  the  external  magnetic  field  Bq^2  [see  Eqs.  (2)  and  (4)  and  Assumption 
(e)].  In  Eqs.  (5)  - (7),  lower  case  jg  denotes  mechanical  momentum. 

Any  distribution  function  that  is  a function  only  of  the  single- 
particle constants  of  the  motion  satisfies  the  steady-state  Vlasov 
equation  (3/3t*0).  For  present  purposes,  we  assume  an  equilibrium 


distribution  function  of  the  form. 


f°(H,pa,p_) 2_2_  fi(u)fi(Pe-P0)  , 


0’  z' 


(8) 


2™Vo 


,2  2. 


where  n^-const.  is  the  electron  density  at  r»Rg,  Pg=-(e/2c) (Rg~a  )Bg= 
const,  is  the  canonical  angular  momentum  of  the  electrons, 


U"H- WYomc  hb 

is  an  effective  energy  variable,  0^=const.  is  defined  by  g^*Vjj/ca 

2 1/2 

(y^-1)  a°d  Yq  and  are  constants.  In  the 

subsequent  analysis,  it  is  shown  that  the  axial  velocity  profile 

associated  with  Eq.  (8)  is  uniform  over  the  beam  cross-section, 


(9) 


with  V (r)-B.  c.  Furthermore,  since  all  of  the  electrons  have  the  same 
z b 

canonical  angular  momentum  Pq=Pq^0,  no  electrons  pass  through  the 
axis  of  symmetry  (r=0) , and  the  equilibrium  electron  density  profile 


0,  , ,3  ,0 

n (r)=  d pf 
e e 


calculated  from  Eq.  (8)  corresponds  to  a hollow  electron 


Several  pertinent  equilibrium  properties  can  be  deduced  for  the 

class  of  thin-beam  equilibria  described  by  Eq.  (8).  For  this  purpose, 

it  is  useful  to  transform  the  energy  variable  U defined  in  Eq.  (9) 

from  momentum  variables  (p  , p0,p  ) appropriate  to  the  laboratory 

r u z 

frame  to  momentum  variables  (p1,  p',  p')  appropriate  to  a frame  of 

r O Z 

reference  moving  with  velocity  B^c^  (the  mean  axial  velocity  of  the 
electron  beam).  Here  ^ is  a unit  vector  in  the  z-direction.  The 

4 

relevant  transformation  is  given  by 


pr«p;,  p0-p;,  Pz=Ybp;+YbY,m6bc,  Y-Yb(Y’+6bP;/mc) 


U ■ ~ (Y'-Y0)mc2  , 0 

yb 

? 2 2 1/2  2 22  1/2 

where  Y“(l+J>  /m  c ) and  y'  = (1+J^'  /me)  . After  some  straight- 

0 3 

forward  algebra  it  follows  from  Eqs.  (8)-(ll)  that  fgd  p transforms 
according  to 


f0,„  p p n0R0(Y,+ebPz/mC) 

fe(H,p0,Pz)d  p-  2i7mY0Y’ 


6(Y’mc2-Y0mc2)6(P’-P0)d3p'  , 


3 3 

where  d p-dp^dp^dp^  and  d p'*dpMp^dp^.  It  is  evident  from  Eq.  (12) 

2 

that  Yginc  -const,  can  be  identified  with  the  total  electron  energy  in 

a frame  of  reference  moving  with  axial  velocity  8bc.  On  the  other 

2 1/2 

hand,  since  the  constant  Yb  is  defined  in  terms  of  $b  by  8b*(Yb~l)  /yb 

2 2 -1/2  2 

[cf.,  discussion  following  Eq.  (9)],  it  follows  that  Ybmc  ”(l_8b)  mc 


is  the  energy  associated  solely  with  the  directed  axial  motion  3^c 

(i.e. , the  energy  calculated  for  pr=p0=pz=O). 

For  a thin-beam  equilibrium  consistent  with  Eq.  (1),  the  energy 
2 

variable  y'mc  in  Eq.  (12)  can  be  approximated  by 

,2.2  2.2-2  2,  ,2.1/2  n 

y me  =c(m  c +m  u>cp  +pL  ) U 


where  use  has  been  made  of  P0=Pq,  p is  defined  by 


P=r-Rn  , 


to  =eB_/mc  is  the  nonrelativistic  electron  cyclotron  frequency,  and 
C U 

2 2 2 

Pi  =p'  +p'  • Substituting  Eq.  (13)  into  Eq.  (12),  and  representing 


Jd3p'=f  dp'[2lTdar  dp]p;  , 

3 J- a.  6Jo  3 0 


it  is  straightforward  to  show  that  the  electron  density  profile  can 
be  expressed  as 


where 


n®(r)=jd3pf°(H,P0,Pz)=no  ®(a2-p2)  , 


t 2 ^V2n,- 

a=(Yn-D  c/u 


is  the  half-thickness  of  the  beam,  and  #(x)  is  the  Heaviside  step 


function  defined  by 


i: : 


x<0  , 


From  Eq.  (16),  we  note  that  the  thin-beam  assumption  (s<<Rq)  is  equivalent 

2 -2  2 2 

to  the  requirement  Yq<<1+(dcRq/c  . Since  the  electron  beam  is  assumed 
to  be  thin,  we  approximate  RQ/r=l  in  Eq.  (15),  and  the  electron  density 
profile  reduces  to 


0,  , (no  ’ lr-Rol<a- 

ne(r)=i 

\ 0 , otherwise 


Making  use  of  v =p  /^“(p'  /m+r ' B.  c)  / (y'tfi,  p'/mc)  , we  obtain  the 

Z Z Z D D Z 

axial  velocity  profile 


„0,  s , ,3  r0w,  ,3  ,0. 
V (r)=(  d pv  f )/(  d pf  ) 
z z e J e 

=8,c=V  =const . , 

D D 


for  |r-RQ|<a.  That  is,  the  axial  velocity  profile  is  uniform  over 
the  beam  cross  section.  In  a similar  manner,  the  azimuthal  velocity 
profile  associated  with  the  equilibrium  distribution  function  in  Eq.  (8) 


can  be  expressed  as 


TT0 . v / ,3  / / ,3  .Ov 

V9(r)-(jd  pve£e)/(jd  pfe) 


YbY0m  ' r 2c  °' 


Making  use  of  Pg=-e (Rg-aZ)Bg/2c  [definition  following  Eq.  (8)]  and  the 
thin-beam  assumption  (a«RQ)  , the  preceding  expression  for  V^(r) 


can  be  approximated  by 


ve(r> 


eBn  a) 

0 c 

p = — p 

Vomc  Yb 


for  |r-Rg|<a. 


Finally,  for  the  equilibrium  distribution  function  in  Eq.  (8), 
the  radial  and  axial  elements  of  the  electron  pressure  tensor  are  defined 


Prr(r)=T°(r)n°(r)=[d3pvrPrf°(H,P0,Pz)  , 


Pzz(r)=Tz(r)ne(r)=jd3p(  W fe(H’P0  ■ Pz>  * 

where  pz=(  d3p  p^  f°)/(  d3p  f °) , vr=pr/ym,  and  T°(r)  and  T^Cr)  are  the 


11 

effective  radial  and  axial  electron  temperature  profiles,  respectively. 

Substituting  Eq.  (12)  into  Eq.  (20),  it  is  straightforward  to  show  that 

a2 

T°(r)=T°(r)=I^-  (a2-p2)  (21) 

for  |r-Rg|<a.  Furthermore,  the  off-diagonal  elements  of  the  equilibrium 
pressure  tensor  are  identically  zero,  i.e.,  Prz(r)=Pzr(r)=0. 

The  self-consistent  electron  density  and  temperature  profiles 
are  illustrated  in  Fig.  2 for  a thin-beam  equilibrium.  Evidently,  for 
the  choice  of  equilibrium  distribution  function  in  Eq.  (8),  the 
electron  density  profile  n^(r)  is  rectangular  with  sharp  radial  boundaries 
[Eq.  (17)],  the  axial  velocity  profile  V^(r)=B,c  is  uniform  over  the 

Z D 

beam  cross  section  [Eq.  (18)],  the  azimuthal  velocity  profile  V^(r)= 
o^p/y^  varies  linearly  over  the  beam  cross  section  [Eq.  (19)], 
and  the  radial  and  axial  temperature  profiles  are  parabolic  [Eq.  (21)]. 
Further  equilibrium  properties  of  intense  hollow  beam  equilibria 
described  by  Eq.  (8)  are  discussed  in  more  detail  in  Ref.  4 . 


Ill . STABILITY  PROPERTIES  FOR  AZIMUTHALLY  SYMMETRIC  PERTURBATIONS 


A.  Linearized  Vlasov-Maxwell  Equations 


In  this  section,  we  make  use  of  the  linearized  Vlasov-Maxwell 
equations  to  investigate  stability  properties  for  azimuthally  symmetric 
perturbations  (3/30*0)  about  the  thin,  hollow-beam  equilibrium  described 
by  Eq.  (8) . We  adopt  a normal-mode  approach  in  which  all  perturbations 
are  assumed  to  vary  with  time  according  to 

<5i|j(x, t)=ij(x) exp{-iojt}  , 


with  lmto>0.  The  Maxwell  equations  for  the  perturbed  electric  and 
magnetic  fields  can  be  expressed  as 


VxE(x)=i  - B(x)  , 

% 'Vf  c % 'V 

(22) 

VxS(x)  =—  J(x)-i  - E(x)  , 

0/  % c 'V/  a/  c ^ V 


where 


(23) 


is  the  perturbed  current  density.  In  Eq.  (23), 

A 

[0  v * xB  (x  * ) 

fe(frZ)Me]^  dT  exp{_iwT}  }'|^  fe  (24) 

is  the  perturbed  distribution  function,  x*t'-t,  and  the  particle 
trajectories  jc'(t')  and  jj'(t')  satisfy  djg'/dt'-v*  and  dp'/dt'  = 

'Xj 

A 

-e^'xBpe^/c,  with  "initial"  conditions  ^c'(t'*t)*jc  and  ^'(t'-t)*^. 

To  make  the  theoretical  analysis  tractable,  we  Fourier  decompose 
the  z-dependence  of  all  perturbed  quantities  according  to 
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= I dkij»k(r)  exp{  ikz  } 


(25) 


where  k is  the  axial  wavenumber.  Note  that  the  perturbation  amplitudes 
in  Eq.  (25)  are  independent  of  9,  because  the  present  stability 
analysis  is  restricted  to  azimuthally  symmetric  perturbations  (2,«0, 
where  J l is  the  azimuthal  harmonic  number).  Making  use  of  Eqs.  (22) 
and  (24) , it  is  straightforward  to  show  that 


h “zk(r>  iW-T-  V'» 


(26) 


where 


2 2.2,2 

p *o)  /c  -k 


and  the  function  <£(r)  is  defined  by 


4>(r)=irE0k(r) 


(27) 


(28) 


Equation  (26)  can  be  expressed  as 

2 

-2-  fi  (r)--SE-S(r)  (29) 

Jr  zr  rui  y ' 

A 

in  the  vacuum  region  outside  the  beam,  where  J , (r)=*0. 

The  perturbed  azimuthal  electric  field  E0k(r)  is  continuous 
across  the  beam  boundaries  (r-R^  and  r~R2)  , as  is  the  function  $(r) 
defined  in  Eq.  (28).  For  azimuthally  symmetric  perturbations, 
it  can  be  shown  that  (Appendix  A) 


|R0(3*/3r)RJ<0[*(R0)]  (30) 

where  0 [<p (Rq)  ] denotes  terms  of  order  $(Rq) . Integrating  Eq.  (26) 
from  r*R^-<S  to  r*R2+6  and  taking  the  limit  6->-0+,  we  obtain  the 
approximate  result 
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Bzk(V->WR 


1 C J»- 


dr  J8k(r) 


(31) 
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where  i|>(R. ) denotes  11m  ±5),  and  use  has  been  made  of  Eqs.  (1) 

J 6-*0  2 

and  ( 30) . 

Similarly,  from  Eq.  (22),  we  obtain 


Making  use  of  the  continuity  of  Ezk(r)  at  r=R^  and  r=R2»  and  the 
inequality  lR0(3E2k/3r)R  |<0[E2k(R0) ] , it  is  straightforward  to  show 
that 


Vek^’-Vek^'T  f-  dr  r j2k(r) 


where  again  use  has  been  made  of  a/Rg«l.  Substituting 


cp 


into  Eq.  (33)  gives 

/ 3E_i 


R2  (ar  ) + ' Rl(ar  ) - = ^ {k[Vrk(R2)_Vrk(Rl)1 

R2  h 


- ip  dr  r J 


c(r)}  . 


Since  the  perturbed  radial  current  density  vanishes  outside  the  beam, 
we  note  that  }rk(V  *^rk^R2^*^'  Therefore  Eq.  (35)  can  be  expressed  as 

r2  - Ri  ( ST*)  - • -1  ^ f - dr  r d,k(r)  • 


For  convenience,  in  the  subsequent  analysis  we  introduce  the 
normalized  electric  and  magnetic  wave  admittances,^  ^ d+  and  b+, 
defined  at  the  inner  and  outer  surfaces  of  the  electron  beam  by 
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and 


d+-trO/Sr)£zk(r)]  /£zk(R2)  , 

R2 

d_-[r(3/3r)Ezk(r)]  , 

R1 

b+fci«k^/tr(8/3r)8.k(r>I«2  • 

b.-Bzk(R‘)/[r(3/3r)Bzk(r)]Ri  . 


(37) 


(38) 


For  a beam  in  vacuum,  Che  values  of  b+  and  d+  depend  In  detail  on  the 
geometric  configuration,  e.g.,  location  of  the  conducting  wall 
(Appendix  B) . It  should  be  noted  from  Eqs.  (22)  and  (25)  that  there 
are  two  main  classes  of  perturbed  fields,  i.e.,  TE  and  TM  perturbations. 
For  the  TE  mode,  the  perturbed  radial  and  axial  magnetic  fields 
(Brk  and  Bzk)  can  be  expressed  as  functions  of  the  perturbed  azimuthal 

A a 

electric  field  EQk.  On  the  other  hand,  for  the  TM  mode,  Efk  and 
B0k  can  be  expressed  in  terms  of  Ezk*  It  is  therefore  convenient  to 
express  the  right-hand  side  of  Eqs.  (31)  and  (36)  as 


^TT  j - dr  50k(r)*xll(u))^(RO)+iROx12(u,)Rzk(RO)  * 
c p 

C-  dr  r dzk^r^*^*2j/w^  • 

R1 


(39) 


where  the  matrix  elements  Xjj 
Making  use  of  Eqs.  (29)-(39) , 
approximate  results 


(oj)  denote  effective  susceptib ilities.^’^ 
it  is  straightforward  to  obtain  the 


(b_-H)+)^(R0)-x11(o))MR0)+iR0x12(<»>)Ezk(R0)  , 
(d.+d+)Ezk(R0)-i[x21(a))/R0]i(R0)+x22(U)Ezk(R0)  . 


(40) 


I 

\ 


! 


' 
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After  some  simple  algebraic  manipulation  that  makes  use  of  Eq.  (40), 
the  dispersion  relation  for  azimuthally  symmetric  perturbations  can  be 
expressed  as 


tb_+b+-x11(u)  ] [d_+d+-X22(u)  J+x12^*21^"°  • (41) 

Evidently,  an  explicit  evaluation  of  the  susceptibilities  X^j (“) 
is  required  for  a detailed  stability  analysis. 


B.  TE  Mode  Dispersion  Relation 

In  this  section,  closed  expressions  are  obtained  for  the  relevant 
electron  susceptibilities  [Eq.  (39)],  and  the  results  are  used  to 
derive  a simplified  TE  mode  dispersion  relation.  We  specialize 
to  the  case  of  a self-consistent  Vlasov  equilibrium  in  which  all 
electrons  have  the  same  canonical  angular  momentum  and  the  same  total 
energy  in  the  moving  frame  [Eqs.  (8)  and  (12)].  To  simplify  the 
right-hand  side  of  Eq.  (24),  use  is  made  of  Eqs.  (6)  and  (9),  and 
the  identities  SU/S^-v-V^^  and  3P0/3pQ-reQ,  where  eQ  is  a unit  vector 
in  the  0-direction.  The  perturbed  distribution  function  can  then  be 
expressed  as 


where 


is  that  portion  of  the  perturbed  distribution  function  that  is  generated 
from  the  TE  perturbation,  and 
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jM  f° 

ek(r’^)“e  au"  j dT  exP^ 1 * -z)-ojt  ] >{  [E 

+izk(r’)(vi“V> 


rk(r,)"ebBek(r,>K 


(43) 


Is  that  portion  of  f ^ that  is  generated  from  the  TM  perturbation. 

As  indicated  in  Eqs.  (42)  and  (43),  the  particle  trajectories, 

r'(x),  9'(t),  z'(t) , in  the  equilibrium  field  configuration  are  required 

in  order  to  evaluate  the  perturbed  distribution  function.  Making  use 

2 2 

of  Eqs.  (5)  and  (6),  and  Pg*Pg“-(e/2c) (R^-a  )Bq  [see  discussion  following 

Eq.  (8)],  it  is  straightforward  to  obtain  the  required  trajectories. 

Assuming  that  the  electron  orbit  passes  through  the  phase-space 

point  (p,p  ) and  (z,p  ) at  time  t'-t,  we  find 
r z 


p '-pcos (u>c/y)  T+(pr/mo)c) sin(o)c/y)  t , 


z'-z+(p/ym)T  , 
z 


(44) 


9’-(o)c/yR0)  [pcos((i)c/y)T+(pr/ma)c)sin(aic/y)T]  , 

where  p^r-R^,  (D^eBg/mc  is  the  nonrelativistic  electron  cyclotron 
frequency  and  use  has  been  made  of  Eq.  (1). 

From  Eqs.  (39)  and  (42),  the  matrix  element  X^(o>)  is 


determined  from 


Xu(.0i<V-T^|*?dr3jk(r)  . 

C P \ 

where  the  perturbed  azimuthal  current  density  corresponding  to  the 
TE  mode  is  defined  by 


(45) 


— ejd3p  v( 


fE 
9 ek 


£“-(r.«) 


(46) 


For  present  purposes  it  is  also  assumed  that 
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|w-kVb-uc/Yb|«u>c/Yb  (47) 

where  wc»eBg/Ygmc  is  the  electron  cyclotron  frequency  in  a frame  of 
reference  moving  with  axial  velocity  Vb-  Within  the  context  of  Eqs. 

(1) , (30),  and  (47)  it  is  valid  to  approximate 


^(r')-ir’E0k  (r')4(R0)  (48) 

and  to  neglect  the  terms  proportional  to  3f||/3P0  in  Eq.  (42),  since  the 

corrections  associated  with  these  terms  are  of  order  y.  Iw-kV.  -oi  /y,  I /to 

b b c b 1 c 

a/Rg  («1),  or  smaller.  Substituting  Eq.  (48)  into  Eq.  (42),  the 
perturbed  distribution  function  in  Eq.  (42)  can  then  be  approximated  by 

-F  / kVh\  3f!  r° 

fek{r’Z)m-±e  V1"'  « ) 3U1  dT  exp{i[k(z:-z)-a)T]}  (49) 

J —00 

where  use  has  been  made  vf  v'*r'0'.  Since  the  variable  U is  an  even 

u 

function  of  pr  [Eq.  (9)],  we  note  from  Eqs.  (46)  and  (49)  that  any 

portion  of  the  time  integral  in  Eq.  (49)  that  is  an  odd  function  of 

Pr  will  automatically  give  zero  when  the  integration  over  pf  is 

carried  out.  Therefore,  we  simply  omit  terms  proportional  to  odd 

“E 

functions  of  p^  when  calculating  fek  from  Eq.  (49).  Making  use  of 
Eqs.  (44)  and  (47),  the  perturbed  distribution  function  for  the 


TE  mode  is  given  by 


e£_ 

2R„ 


V(R0} 

YOJ-u  -kp  /m 
c z 


(50) 


A 

where  oj^eBg/mc. 

f 3 -E 

The  integral  d p v0  fgk  is  required  to  evaluate  the  perturbed 
azimuthal  current  density.  For  a thin,  hollow  beam,  we  note  from 
Eq.  (44)  that  v *r§  can  be  approximated  by 


vQ*  pwc/y 


(51) 
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Substituting  Eqs.  (8),  (50),  and  (51)  Into  Eq.  (46),  making  use  of 
the  transformations  In  Eq.  (10),  and  carrying  out  the  required 
integration  over  momentum,  the  perturbed  azimuthal  current  density 


corresponding  to  the  TE  mode  can  be  expressed  as 
2 2 


-E  e n0p 

Jek(r)  ' 2Vo““o 


(l  - AfR  I 6(p-a)+6(p+a) 

V u>  ) *(V  j a(a)-kVb-a)^/Yb) 


u^ct((i)-kVb)  (u-kVb-a)c/yb)-(a)eb-kc)2]«(a2-p2) 

" r 2,  , T >2  . , ,2  2.  2 2.  . 3/2 

[c  (a)-kvb-u)c/Yb)  -(o)6b-kc)  u>c(a  -p  )] 

where  is  the  electron  cyclotron  frequency  in  the  moving 

frame.  Substituting  Eq.  (52)  into  Eq.  (45),  we  find 


(52) 


xu<*) 


wj?  Cu-kv 
_E  b 


Yb  c2p2(kc-wBb)2R0 


x£n 


Cft  +■  b) 

C 

cO  - a) 

c 


c[n-(kc-a>Bb)  /(u-kvb)] 

2lo  (kc-o>6v) 

C D 


a(kc  - o)Bb) 


a(kc  - o)8b) 


- 


(53) 


2 2 

where  u)p"4irnQe  /Yq®  is  the  electron  plasma  frequency  squared  in  the  moving 
frame,  and  fi  is  the  Doppler  shifted  frequency  defined  by 


0*w-kV,-u>  / y v • (54) 

The  following  point  of  verification  is  noteworthy.  If  instead 
we  make  use  of  the  equilibrium  distribution  function  (characterized 
by  V°-0) 

f>-V  ‘ «<“-vc2wW  • 

and  then  make  the  appropriate  relativistic  transformation  to  a 
frame  of  reference  moving  with  axial  velocity  “®bc^z»  then  the 
expression  obtained  for  x^tui)  la  identical  to  Eq.  (53). 

Following  similar  algebraic  manipulations,  we  have  evaluated 
the  remaining  susceptibility  matrix  elements  x^’  X21*  and  X 22 ' 


However,  after  some  straightforward  but  tedious  algebra,  it  is  found 
that 


4 I « 

1 

X21 

x22 

Q 

. 

I 

X11 

9 

X11 

u 

ce 

(55) 


Therefore,  within  the  context  of  the  inequalities  in  Eq.  (55),  the 


TE  mode  dispersion  relation  obtained  from  Eq.  (41)  can  be  approximated 
by 

b_+b+a‘X22^u)^  * (56) 

In  Sec.  IV,  we  make  use  of  the  approximate  dispersion  relation  in 
Eq.  (56)  to  investigate  TE  mode  stability  properties  for  azimuthally 
symmetric  perturbations  about  the  hollow-beam  equilibrium  described 
by  Eq.  (8) . 
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IV.  TE  MODE  STABILITY  ANALYSIS  FOR  AZIMUTHALLY  SYMMETRIC  PERTURBATIONS 


Within  the  context  of  Eq.  (55),  the  TE  mode  is  the  dominant 
unstable  perturbation.  In  this  section,  we  make  use  of  Eqs.  (53) 
and  (56)  to  investigate  TE  mode  stability  properties  in  parameter 


regimes  of  experimental  interest 


21-24 


The  number  of  electrons  per  unit  axial  length  can  be  expressed  as 
Ne*4irnQaRg.  Eliminating  n^  in  favor  of  Ng  in  Eq.  (53)  and  substituting 
Eq.  (53)  into  Eq.  (56),  the  TE  mode  dispersion  relation  is  given  by 


v (aHcV 

b +b  = — j 

YbY0  p2Rg(kc-u>Bb)2 


f2+Bg(kc-a)Sb) 


[fi-(kc-m6b)  /(co-kVb)] 
26q  (kc— co6b) 


(fi-B0(kc-u>eb)jj 

22  2 1/2 
where  v=N&e  /me  is  Budker's  parameter,  Bq=(1-1/yq)  , and  use  has  been 

made  of  Eq.  (16).  The  growth  rate  uK=Imai  and  real  oscillation 

frequency  m^Rern  can  be  determined  from  Eq.  (57)  for  a broad  range 

of  system  parameter,  by  solving  numerically  the  full  transcendental 

dispersion  relation.  We  reiterate  that  Eq.  (57)  is  valid  only  when 

I n | = | aj-kV, -a)  /y,  | < <u>  /y,  [Eq.  (47)].  Moreover,  it  should  be  noted 
b c b c b 

that  the  magnetic  wave  admittances  b+  defined  in  Eq.  (B.3)  are  generally 
complicated  functions  of  the  eigenf requency  ui-w^+ia)^. 
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A.  Limit  of  a Tenuous,  Infinitesimally  Thin  Beam 


To  simplify  the  expressions  for  the  magnetic  wave  admittances 
b+  defined  in  Eq.(B.3),  it  is  useful  to  consider  the  limit  where  the 
electron  beam  is  infinitesimally  thin. 


a/R0->°  » 


and  the  beam  density  is  very  tenuous,  i.e. , 


v/Yq-0  . 


Making  use  of  Eq.  (58),  it  is  straightforward  to  show  from  Eq.  (B.3) 
in  Appendix  B that  b_+b+  can  be  approximated  by 

2J1(pRc)Mp2Rj 

b~+V"  J1(pR0)[J1(pR0)N1(pRc)-J1(pRc)N1(pR0)]  » (60) 

2 2,22 

where  p =u>  /c  -k  , J^(x)  and  N^(x)  are  Bessel  functions  of  first  and  second 
kind,  respectively,  and  use  has  been  made  of  N^(x) J^(x)-Jq(x)N^(x)=2/ttx. 

In  the  absence  of  a beam,  the  vacuum  TE  mode  dispersion  relation  is 


given  by 


2^2  ,2,2,  2,_2  2 

p Rc=(m  /c  -k  )Rc=«0n  , 


where  ct-  is  the  nth  root  of  J,(a„  )=0.  For  a very  tenuous  beam 
un  ± un 

[Eq.  (59)],  it  is  evident  from  Eqs.  (57)  and  (60)  that  Eq.  (61)  is 

a good  approximation  to  the  dispersion  relation  in  Eq.  (57).  [The 

right-hand  side  of  Eq.  (57)  describes  the  beam-produced  modifications 

to  the  vacuum  dispersion  relation  (61).]  Taylor  expanding  Eq.  (60) 

2 2 2 2 2 

about  (m  /c  -k  )R  *a_  , it  is  straightforward  to  show  that  b +b,  can 
c Un  - + 

be  approximated  in  leading  order  by 

1 IW  ( »2  .2  jW\  fJ2(a0n)  } . 

- + ' 2p2  1R0  ) L2  ' ‘ R \>  ' 
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i 


Substituting  Eq.  (62)  into  Eq.  (57),  we  obtain  the  approximate 


dispersion  relation 
a)2  . 2 a0n 

~2~  k ~~T 

c R 

c 


2v 

Vo 


J1(«n_Rn/R Jl2  /m-kVu  \2 


1 On  0 c 


RcJ2(a0n) 


n-(kc-u)$b)  / (w-kV^) 
1 260(kc-a>$b) 


kc-wB, 


0+SQ(kc-w6b) 


(63) 


n-S0(kc-o)6b) 


for  a tenuous,  infinitesimally  thin  beam. 

To  lowest  order,  the  eigenfrequency  u>  and  axial  wavenumber  k are 
obtained  from  the  simultaneous  solution  of  Eq.  (61)  (the  vacuum  TE 
mode  dispersion  relation)  and  the  condition  for  cyclotron  resonance 
[Eq.  (47)] 


<i>=kV,-ko  ly,  . 

D C D 


(64) 


In  addition,  as  shown  in  Appendix  C,  the  region  of  k space  corresponding 
to  instability  is  very  narrow-band  with 


kc=ujf3. 


b * 


(65) 


Solving  Eqs.  (64)  and  (65)  for  the  characteristic  frequency  and  wavenumber 
(a),k)=(a>0,k0) , we  find  (Fig.  3) 


Wb  ’ 


Wb6b/C  * 


(66) 


2 -1/2 

where  use  has  been  made  of  Yb=(l-&b)  . Since,  for  a very  tenuous 

beam,  we  also  require  that  (u)g,kg)  solve  Eq.  (61)  in  leading  order, 
it  follows  that  instability  exists  only  in  the  narrow  region  of 
parameter  space  satisfying 


i 
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1 


On 


(67) 


This  feature  is  also  evident  from  the  numerical  studies  of  the  dispersion 
relation  presented  in  Sec.  IV. B.  Finally,  for  future  reference, 
it  should  be  noted  from  Eqs.  (61)  and  (65)  that  the  group  velocity 
can  be  proximated  by 

2 


v = dm  = kc_  a v 
g dk  a)  b 


(68) 


That  is,  for  a very  tenuous  beam,  the  group  velocity  of  the  unstable 
waves  is  approximately  equal  to  the  beam  velocity. 

Appendix  C and  the  previous  paragraph  give  some  useful  insight 
into  the  region  of  parameter  space  corresponding  to  instability. 

Most  important  is  the  fact  that  the  instability  is  narrow-band 


in  k-space  with  | kc-ooB^ | <<c/Rq  [Appendix  C and  Eq.  (65)].  Expanding 


the  logarithmic  contribution  for  small  |kc-ajB^|,  the  dispersion  relation 


in  Eq.  (63)  can  be  expressed  in  the  approximate  form 


2 9 <*n 

a)  ^2  _ On 


2v 


J, («n-*n/R  ) 


1 On  0 c; 


YbY0  *■  RcJ2(a0n) 


rio-kV, 


[o)-kVb-a,c/Yb 


602(u-kVb)2 


(69) 


3(u-kVb-(jjc/Yb)' 


2 1/2 

where  0q=(1-1/Yq)  • The  general  features  of  Eq.  (69)  are  qualitatively 
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similar  to  the  result  obtained  by  Friedman  et  al.  for  perturbations 


about  a self-consistent  Vlasov  equilibrium,  f ^=const . xg (H^Ygmc2)  x 


6 ( p -Y.  mg,  c) 6(P  -Prt) , where  all  electrons  have  the  same  perpendicular 

Z D D “ U 


energy  H^-YqJHC  . As  indicated  in  the  introduction,  however,  Eq. 
(69)  differs  in  detail  as  do  the  corresponding  stability  results. 


For  example,  the  last  term  on  the  right-hand  side  of  Eq.  (69)  is 


2 2 2 2 2 2 

proportional  to  8g(o)-kVb)  /3,  instead  of  Bq(o)  -k  c )/2  as  obtained  in 


i 
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Ref.  21.  Therefore,  depending  on  the  region  of  parameter  space,  the 

stability  properties  obtained  from  Eq.  (69)  can  differ  considerably 

from  the  results  obtained  in  previous  analyses  of  the  cyclotron  maser 
21 

instability.  Moreover,  the  dispersion  relation  in  Eq.  (63)  predicts 

instability  for  a very  narrow  range  of  k-values , whereas  previous 
21 

studies  have  not  been  able  to  identify  the  range  of  unstable  k values 
in  a consistent  manner. 

Finally,  the  necessary  and  sufficient  condition  for  Eq.  (69)  to 
have  unstable  solutions  (lmw>0)  is  given  by 


„4  4 

60  > 3 


VC 


[Jl<“onR0/V] 


2 

Y0Yb“c 


RcJ2(o0n> 


(70) 


which  is  qualitatively  similar  to  the  result  obtained  by  Sprangle 
and  Drobot^  for  a non-self-consistent  equilibrium  slab  configuration. 


B.  Numerical  Analysis  of  Dispersion  Relation 


In  this  section  we  make  use  of  the  dispersion  relation  (57) 

to  carry  out  a numerical  investigation  of  stability  properties.  In 

this  regard,  no  a priori  assumption  is  made  that  the  beam  is  very  tenuous 

(as  in  Sec.  IV. A) . However,  use  is  made  of  the  fact  that  the  narrow- 

band  of  unstable  k values  satisfies  | kc-wS^ I <<c/Rq.  Expanding  the 

2 

logarithm  in  Eq.  (57),  and  evaluating  k,  p , and  b++b_  at  (u),k)  = 

(oi0,k0)  [Eq.  (66)],  the  dispersion  relation  (57)  can  be  approximated  by 


p (b_+b+)+(io-a>0) 


(b-+v 


(m-Vb) 


, 2 (w-k  V, -w  /y  ) 

w0»ko  YbY0R0  0 b c b 


, . fo  (“-koV 
■ 


(71) 


3 (u,-k0Vb-u,c/yb) 

*2  2 2 2 a a 

where  p Hu)„/c  -k_,  and  b_+b  =[b_+b  ] , 

U U — -r  — , Kq 


is  defined  by  Eq.  (B.3). 
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The  growth  rate  ui^Imw  and  real  oscillation  frequency  cor=Reu)  have 

bee*,  calculated  numerically  from  Eq.  (71)  for  a broad  range  of  system 

parameters  a/RA,  R /Rrt,  v,  6,  and  For  R /R.  corresponding  to 

UcU  d u c U 

A A 

b_+b+“0,  we  have  Taylor  expanded  (b_+b+)  about  co=coq , approximating 

(b  +b.)  in  leading  order  by  (b  +b , ) = (u)-a)A)  [d(b  +b^)/da)]  , . For 

- + - + 0 - + <*>0>k0 

present  purposes,  to  illustrate  the  influence  of  the  geometric  configuration 
on  stability  behavior,  we  calculate  the  normalized  magnetic  wave  admittances 
b+.  Shown  in  Fig.  4(a)  is  a plot  of  b++b_  versus  rc/rq*  for  Bq=0. 2, 

6^=0.143  and  a/RQ=0.1.  It  is  evident  from  Fig.  4(a)  that  b_+b+*0  for 
Rc/Rq“2.18  and  Rc/Rq=3.75,  which  corresponds  to  the  TE  mode  dispersion 
relation  without  beam-produced  modifications.  Figure  4(b)  shows  a plot 
of  the  normalized  growth  rate  versus  Rc/Rq  obtained  from  Eq.  (71) 

for  v*0.001  and  parameters  otherwise  identical  to  Fig.  4(a).  Two 
important  features  are  noteworthy  from  Fig.  4.  First,  for  a low  beam  den- 
sity (v<0.001),  the  maximum  growth  rate  occurs  at  the  value  of  R /R„ 

'v  c 0 

for  which  b_+b+  is  equal  to  zero.  For  example,  in  Fig.  4(b),  the 
system  is  most  unstable  at  Rc/Rq*2.18.  Second,  for  a low  beam  den- 

A A 

sity,  the  system  is  unstable  only  in  the  range  where  b +b+>0  [compare 
Figs.  4(a)  and  4(b)]. 

Shown  in  Fig.  5(a)  is  a plot  of  (b  +b+)  versus  rc/Rq  for  Bq=0. 4, 
6^*0.286  and  a/Rg=0.2.  Also  shown  in  Fig.  5(b)  are  plots  of  the 
normalized  Doppler  shifted  real  frequency  Re(J2/o)c)  = (a)r-kV^-ajc/Y^)  /a>c 
and  growth  rate  u^/a)(  versus  rc/rq  ^or  v=0.001  and  parameters  otherwise 
identical  to  Fig.  5(a).  Note  that  Re(ft/coc)  is  plotted  only  for  the 
range  of  Rc/Rg  corresponding  to  instability.  Evidently,  for  the  low 
beam  density  in  Fig.  5,  the  stability  properties  are  qualitatively 
similar  to  those  in  Fig.  4.  However,  it  is  also  evident  from 
Figs.  4(b)  and  5(b)  that  the  growth  rate  increases  rapidly  as  8q  is 
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increased.  Moreover,  we  note  from  Fig.  5(b)  that  Refi/wc  assumes 

a maximum  value  (0.019)  at  Rc/Rq=2.375  and  decreases  rapidly  to  zero 

when  R /Rn  is  increased, 
c u 

We  conclude  this  section  by  emphasizing  that  stability  properties 

also  exhibit  a sensitive  dependence  on  beam  density.  This  is  illustrated 

in  Fig.  6 where  the  normalized  growth  rate  is  plotted  versus 

Rc/Rq  for  two  different  values  of  v,  and  parameters  otherwise  identical 

to  Fig.  5(a).  For  intermediate  or  high  beam  densities  (v>0.01), 

the  maximum  growth  rate  no  longer  coincides  with  values  of  R^/Rq 

for  which  (b_+b+)=0.  For  example,  in  Fig.  6,  the  maximum  growth  rate 

occurs  for  R /RA=2.45  when  v=0.01,  and  for  R /R  =2.51  when  v=0.1. 
c 0 c c 

Evidently,  the  value  of  Rc/Rg  corresponding  to  maximum  growth  rate 
is  a slowly  increasing  function  of  beam  density.  However,  the  value 
of  the  maximum  growth  rate  increases  relatively  rapidly  with 
increasing  beam  density.  We  also  note,  for  high  beam  densities, 

A /V 

that  the  system  can  be  unstable  even  for  negative  values  of  (b  +b+) 

(Fig.  6).  This  is  not  the  case  for  low  beam  densities  (Fig.  4). 
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V.  CONCLUSIONS 


In  this  paper  we  have  formulated  a self-consistent  theory  of  the 

cyclotron  maser  instability.  The  stability  analysis  was  carried  out 

for  an  infinitely  long  hollow  beam  propagating  parallel  to  a uniform 

magnetic  field  bq£z»  within  the  framework  of  the  linearized  Vlasov- 

Maxwell  equations.  The  equilibrium  configuration  and  basic  assumptions 

were  summarized  in  Sec.  II. A.  In  Sec.  II. B,  equilibrium  properties 

were  calculated  for  the  choice  of  electron  distribution  function  in 

which  all  electrons  have  the  same  value  of  canonical  angular  momentum 

2 

(Pq)  and  the  same  value  of  energy  (Yginc  ) in  a frame  of  reference 
moving  with  axial  velocity  8,c  [Eq.  (8)].  Stability  properties  for 

D 

azimuthally  symmetric  perturbations  were  examined  in  Secs.  Ill  and  IV. 
As  shown  in  Sec.  III.B,  the  TE  mode  corresponds  to  the  most  unstable 
perturbation.  Moreover,  the  analysis  in  Appendix  C predicted  that 
instability  occurs  for  a very  narrow  range  of  k values  satisfying 
|kc-8^u)|<<c/RQ.  A detailed  numerical  analysis  of  the  approximate 
dispersion  relation  (71)  was  presented  in  Sec.  IV. B.  One  of  the  most 
important  conclusions  of  this  study  is  that  the  geometric  configuration 

A A 

(b_+b+)  has  a large  influence  on  stability  behavior,  particularly 
when  the  beam  density  is  low.  Moreover,  the  maximum  growth  rate 
increases  relatively  rapidly  with  increasing  beam  density. 
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APPENDIX  A 


EVALUATION  OF  PERTURBED  AZIMUTHAL  ELECTRIC  FIELD 


From  Eqs.  (22)  and  (28),  we  find  that 


d * 


- f 


(A.  1) 


Since  the  perturbed  azimuthal  current  density  is  even  function  of  p 
[Eq.  (52)],  it  is  straightforward  to  show  that 


Vv-t  i“2t(4>+52k<vj  • 


(A.  2) 


Making  use  of  Eqs.  (29),  (38),  (A.l)  and  (A. 2),  we  obtain 


RQ (d^/dr) Rq/$(R0)  - y P2RQ(b+-b_) 


(A.  3) 


From  Eq . (A. 3) , it  follows  that 


|R0(d*/dr)RJ<0[$(R0)]  , 


provided  |p  RQ(b+-b_) |<1. 


(A. 4) 
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APPENDIX  B 

EVALUATION  OF  MAGNETIC  WAVE  ADMITTANCES 

In  this  section,  we  obtain  expressions  for  the  wave  admittances 
at  the  boundaries  of  a hollow  electron  beam  in  a cylindrical  waveguide. 

Since  the  perturbed  current  density  vanishes  in  the  vacuum  region 
outisde  the  beam.  Maxwell's  equation  in  this  region  can  be  expressed 
as 

' h 1 h + <*•» 

2 2 2 2 

where  p *u)  /c  -k  . In  obtaining  Eq.  (B.l),  use  has  been  made  of  Eqs. 

(29)  and  (A.l).  The  solution  to  Eq.  (B.l)  is  given  by 

{ AJ  (pr)+BN_(pr)  , R.<r<R  , 

B .(r)  - < 0 0 1 C (B.2) 

* CJQ(pr)  , 0<r<R1  , 

where  J^(pr)  and  N^(pr)  are  Bessel  functions  of  the  first  and  second 

kind,  respectively,  and  R is  the  radius  of  the  waveguide  (Fig.  1). 

c 

In  Eq.  (B.2),  the  constants  A,  B,  and  C are  related  by  Eq.  (28)  and 

the  boundary  condition  [(3ldr)B^(r)  ]R  =0.  After  some  straightforward 

c 

algebra,  the  magnetic  wave  admittances  defined  in  Eq.  (38)  can  be 
expressed  as 

! J0(pR2)N1(pRc)-J1(pRc)N0(pR2) 

b+  * R^p  J1(pR2)N1(pRc)-J1(pRc)N1(pR2)  » 

(B.  3) 

. 1 J0(pRl) 

- RlP  J1(pR1)  * 

where  R^  and  R2  are  the  inner  and  outer  radii  of  the  beam,  and  use  has  been 
made  of  dJg(x) /dx=-J^(x)  and  dNg(x) /dx=-N^(x) . Closed  expressions  for  the 
electric  wave  admittances  d+  [Eq.  (37)]  can  be  obtained  in  a similar  manner' 
but  are  not  required  in  the  TE  mode  stability  analysis  in  Sec.  IV. 
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APPENDIX  C 

NYQUIST  STABILITY  ANALYSIS  FOR  PERTURBATIONS  WITH  k^B,  u>/c 

D 

In  this  section,  we  make  use  of  the  Nyquist  method  to  obtain 

necessary  and  sufficient  conditions  for  instability  from  the  TE  mode 

15  26 

dispersion  relation  [Eq.  (57)].  * The  analysis  is  carried  out 

for  perturbations  with  the  axial  wavenumber  k satisfying 

q“k-Bbu)/c^0  . (C.l) 

I 

Consistent  with  Eq.  (47),  for  present  purposes  we  approximate  w~kVb+u)c/Yb 

2 2 2 2 

in  the  definitions  of  q [Eq.  (C.l)],  p *u>  /c  -k  , and  the  magnetic  wave 
admittances  b+  [Eq.  (B . 3) ] . Assuming  (b  +b+)^0,  and  defining 

z-0/B0qc»(to-kVl5+ioc/Yb)/BQqc  , 

s-qcYb/B0o)c  , (C. 2) 

2 22  3,  2 

a-c  p R0Y0Yb/vuic  » 

it  is  straightforward  to  express  the  TE  mode  dispersion  relation  in 
Eq.  (57)  as 

D(z)-(z-s)£n(|^)  - 2 * -2(b_+b+)q2a  , (C.3) 

where  the  complex  variable  z is  defined  in  Eq.  (C.2).  Note  from 
Eq.  (C.3)  that  the  dispersion  relation  is  an  even  function  of  q so  that 
stability  properties  do  not  depend  on  the  sign  of  q.  Therefore,  we 
assume  q>0  without  loss  of  generality. 

Nyquist  plots  of  the  function  D(z)  defined  in  Eq.  (C.3),  which 
map  the  upper  half  z plane  (contour  f)  onto  the  complex  D plane 
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(contour  I”),  enable  us  to  study  the  stability  properties  predicted  by 

Eq.  (C.3)  [Figs.  7 and  8].  Figures  7(a)  and  8(a)  are  plots  of  contour 

T for  s<l  and  s>l,  respectively.  Branch  cuts  originating  from  the 

branch  points  at  z*±l  must  be  introduced,  as  shown  in  Figs.  7(a)  and 

8(a),  to  make  the  function  D(z)  single  valued.  Plots  of  D(z)  are 

shown  in  Fig.  7(b)  for  s<l,  and  in  Fig.  8(b)  for  s>l.  In  obtaining 

the  contour  T'  in  Figs.  7(b)  and  8(b),  use  has  been  made  of  the 

definition  £nz“£nr+i0  (— tt<9 <tt)  , where  z=rexp{10}. 

The  Nyquist  analysis  indicates  that  the  system  is  stable  (lmft<0) 

provided  the  contour  T'  in  the  complex  D plane  does  not  enclose  the 
2 

point  -2(b_+b+)q  ot.  Therefore,  it  follows  from  Figs.  7(b)  and  8(b) 
that  a necessary  condition  for  instability  is 

s <1  . (C.4) 

| 

Moreover,  from  Fig.  7(b),  we  also  note  that 

0<(b_+b+)q2a<l  (C.5) 

2 

is  a necessary  condition  for  instability.  Since  q a>0,  we  conclude 
that  (b  +b+) >0  is  required  for  instability,  which  is  consistent 
with  the  results  obtained  from  a numerical  analysis  of  the  dispersion 
relation  (see  Fig.  4).  Combining  the  two  inequalities  in  Eqs.  (C.4) 
and  (C.5),  and  making  use  of  Eq.  (C.2),  it  is  straightforward  to  show 
for  typical  beam  parameters  that 

| q | ■ | k-8^a>/ c | «1/Rq  (C.6) 

is  required  for  instability,  which  corresponds  to  a very  narrow 
bandwidth  in  k-space.  The  inequality  in  Eq.  (C.6)  can  also  be 
demonstrated  for  the  case  (b_+b+)=0,  by  Taylor-expanding  (b_+b+) 
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] 

about  w-kV^-H^/Y^*  keeping  first-order  terms  in  ft,  and  making  use  of 
[d(b_+b+)/do>] ky  / >0. 

D C D 
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FIGURE  CAPTIONS 

Fig.  1 Equilibrium  configuration  and  coordinate  system. 

Fig.  2 Electron  density  [Eq.  (17)]  and  temperature  profiles  [Eq.  (21)]. 

Fig.  3 The  straight  lines  w=kV  -ho  / y,  and  w=kc/6,  intersect  at  (u)n,k,»)  = 

D C D D U U 

2 2 2 2 2 1/2 

Yt/c).  The  curve  a)=(c  k +a_  c /R  ) passes 
CD  D C D CJn  c 

through  (wgjkg)  provided  a0nC^Rc=u,c ' 

Fig.  4 (a)  Plot  of  sum  of  magnetic  wave  admittances  (b  +b+)  (Appendix  B) 

versus  R /R0,  for  0 =0.2,  g =0.143  and  a/R  =0.1.  (b)  Plot  of 

normalized  growth  rate  versus  Rc/RQ  [Eq.  (68)],  for  v=0.001 

and  parameters  otherwise  identical  to  Fig.  4(a). 

Fig.  5 (a)  Plot  of  sum  of  magnetic  wave  admittances  (b  +b+)  (Appendix  B) 

versus  Rc/Rq,  for  8Q=0.4,  gb=0.286  and  a/RQ=0.2.  (b)  Plot  of 

the  normalized  Doppler  shifted  real  frequency  Ren/mc=Re(iD-kVb-uc/Yb) /u>c 
and  growth  rate  u^/u^  versus  Rc/Rq  [Eq.  (68)],  for  v-0.001  and 
parameters  otherwise  identical  to  Fig.  5(a). 

Fig.  6 Plot  of  normalized  growth  rate  u>j/u>c  versus  Rc/Rq  [Eq.  (68)], 
for  v=0.01  and  v=0.1,  and  parameters  otherwise  identical  to 
Fig.  5(a). 

Fig.  7 Map  of  upper-half  z plane  (contour  T)  onto  the  D plane  (contour  r') 
for  s<l. 

Fig.  8 Map  of  upper-half  z plane  (contour  f)  onto  the  D plane  (contour  T') 


for  s>l. 


